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ABSTRACT 



This thesis examines radio propagation conditions over trans-equatorial (TE) paths. 
The study precedes Project PENEX, a field experiment to measure and collect calibrated 
HF skywave signal strength data for polar, equatorial, and near-vertical incidence 
propagation paths. PENEX will benchmark the absolute accuracy of the signal-to-noise 
models in the MEDUSA propagation model now being developed by the Naval 
Command, Control and Ocean Surveillance Center. 

Since only minimal information is available about TE paths, a comprehensive 
review of the published literature on TE experiments was completed, with emphasis on 
TE paths between magnetic conjugate points. Data from such paths have revealed the 
presence of unusual propagation modes which are not predicted by standard propagation 
programs such as PROPHET, IONCAP, and AMBCOM. The review of the literature 
revealed that Stanford Research Institute (SRI) published measured data on a TE path 
between the Pacific islands of Kauai and Rarotonga. 

A comparison was made between the SRI data and predictions for the same path 
to assess the usefulness of current prediction programs for TE paths. The SRI AMBCOM 
program was used for this comparison. As expected, sizeable differences were found 
between the predicted and measured results, especially during times when unusual 
propagation modes were present. This suggests that prediction programs should be 
modified to include the observed TE modes. 
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I. INTRODUCTION 



A. GENERAL 

Radio propagation aspects of trans-equatorial (TE) paths are examined in this thesis. 
This study precedes a series of actual propagation measurements under Project PENEX 
to be conducted in 1993-94 by the Naval Security Group Command, the Naval 
Command, Control, and Ocean Surveillance Center, and the Naval Postgraduate School. 

Since very little information is available in textbooks and standard documents about 
HF TE propagation, a comprehensive review of the literature was completed. A 
bibliography of pertinent documents is provided at the end of the thesis. Information 
obtained from documents in the bibliography were of interest and value in identifying TE 
information and experiments. 

A number of experiments have been conducted over trans-equatorial paths, and 
many of the researchers involved have reported unusual TE propagation conditions. 
These unusual conditions have included propagation of unusual modes, at unexpected 
frequencies and anomalously high signal amplitudes. One particular experiment was of 
special interest, since one of the paths involved will most likely be the TE path for the 
PENEX study. Stanford Research Institute (SRI) conducted a series of experiments on 
TE paths in the Pacific region during the nuclear test series of 1962. The path from 
Kauai in the Hawaiian Islands to Rarotonga in the Cook Islands was of special interest 
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because the two islands are situated very nearly at magnetic conjugate locations. Sounder 
data from this path were published by SRI [Ref. 1 ]. 

A comparison was made between this sounder data and a standard propagation 
prediction program. While several prediction programs could have been used for this 
purpose (e.g.. PROPHET, IONCAP and others), the SRI AMBCOM program was used 
in this study. 

This thesis is divided into five chapters. The first chapter provides an overview of 
Project PENEX, and a brief description of the morphology of the TE ionosphere as it 
applies to the AMBCOM program. The second chapter will describe the AMBCOM 
algorithms. Chapter III will characterize the 1962 SRI measured data and the 
experimental setup. The fourth chapter will compare AMBCOM predictions with SRI's 
measured data. Conclusions and recommendations are in Chapter V. The interested 
researcher will find a bibliography of the equatorial ionosphere and TE propagation at the 
conclusion. 

B. PROJECT PENEX 

PENEX is an acronym for Polar. Equatorial, and NVIS (Near-Vertical Incidence 
Skywave) Experiments. The objective of the PENEX experimental program is to measure 
and collect calibrated HF skywave signal strength data for the purpose of benchmarking 
the absolute accuracy of the signal to noise (S/N) models in the MEDUSA propagation 
models, now being developed by the Naval Command, Control, and Ocean Surveillance 
Center. Based on ten years operational experience and numerous comparisons to 
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experimental data, the radiowave propagation models in PROPHET and MEDUSA 
provide very accurate median predictions for mid-latitudes. The same propagation models 
perform only marginally in predicting short term signal variations. In fact none of the 
existing propagation codes in use do a good job. 

PENEX researchers will target the propagation characteristics in auroral and polar 
cap regions, trans-equatorial (TE) regions, and the propagation mode known as Near- 
Vertical Incidence Skywave (NVIS). Project PENEX is a 2-year study (1993-1994) of 
HF propagation, unique in utilizing a wide-band spread spectrum matched filter technique. 
Employing direct-sequence spread spectrum modulation techniques, GPS location data, 
and a rubidium clock, researchers expect the data generated to be very highly correlated. 
The wideband signal (about 40 kHz) approach is attractive and offers the following 
features: 

1. Absolute signal recognition in almost any kind of interference environment (spread 
spectrum processing gains > 800). 

2. Sufficient time resolution (12.5 ps time delay resolution) to identify each mode of 
propagation and the power density in that mode. 

3. Significantly reduced output power requirements (< 100W). 

The TE portion of the experiments, the subject of the current thesis, is scheduled 
to be conducted in 1993-94 with a transmitter located on Kauai in the Hawaiian Islands 
and a receiver on Rarotonga in the Cook Islands. 
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c. 



THE TRANS-EQUATORIAL IONOSPHERE 



1. Geomagnetic and Solar Parameters 

a. The Geomagnetic Coordinate System 

The Earth's magnetic field may be approximated by an earth-centered 
dipole directed southward and inclined at about 11.5° to the earth's rotational axis. 
Presently, the northern pole of the dipole is located approximately at 81° N, 84.7° W 
using the geographic coordinate system [Ref 2: p.60] When studying or modeling the 
ionosphere, the geomagnetic coordinate system is commonly used to map the ionosphere 
to an earth-bound coordinate system. The geomagnetic coordinate system is based upon 
the geographic location of the earth's magnetic poles. The longitudinal origin for this 
system is the meridian line which passes through the north and south geomagnetic poles 
and through the geographic south pole. The two systems are related with the equations 
[Ref 3: p.40] 

sin <E> = sin<[> sin(f> 0 + cos<t> cos<}) 0 cos(A. - X 0 ), (1) 

and sin A = cost}) sin(A. - X 0 )/ cos 0, (2) 

where 

4> 0 = geographical latitude for the northern geomagnetic pole, 

X 0 = geographical longitude for the northern geomagnetic pole, 

4> = geographical latitude, 

X - geographical longitude, 

O = geomagnetic latitude and 
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A = geomagnetic longitude. 

Although the AMBCOM program makes the geographic-to-geomagnetic coordinate 
conversion for the user, an understanding of the geomagnetic coordinate system is useful 
in understanding the model. 

b. Geomagnetic Activity Indices 

AMBCOM uses the three-hour Kp index of worldwide magnetic 
disturbance to specify the current state of the earth's magnetic field. This index is based 
upon local K indices which are quasi-logarithmic values prepared at twelve selected 
observatories worldwide to describe the condition of the planetary magnetic field at each 
site. These local K values are corrected to calculate the planetary Kp index. The K r 
index is calculated at three hour intervals for each of eight periods per day. The Kp 
indices range from zero, the least disturbed state, to nine which represents the most 
disturbed magnetic field. The Kp indices used for this thesis are published in the Journal 
of Geophysical Research [Refs. 4, 5, 6, 7] and are listed in Appendix B. For this study the 
indices ranged from zero to seven, indicating that the geomagnetic field varied from a 
nondisturbed state to a considerably disturbed state. 

a Sunspot Number 

Sunspots, related to the solar flux, are characterized by strong magnetic 
fields which may approach 0.4 Tesla. They have approximately an 11-year periodicity 
in occurrence. Their occurrence was measured by the Wolf, or Zurich sunspot number, 
R 2 [Ref. 2: p.29]. 
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R, = k(10g + s) 



( 3 ) 



where 

g = number of sunspot groups observed, 
s = number of observed individual spots and 
k = correction factor. 

Daily values of R z used in the present study were also derived from the 
Journal of Geophysical Research [Refs. 4, 5, 6, 7] and are listed in Appendix B. The R z 
number was discontinued in 1981 in favor of the International Sunspot Number, Rj [Ref. 
2: p. 44], 

2. Ionospheric Layers 

Radio waves may be refracted or reflected as they encounter ionospheric layers 
during propagation. As rays move through various areas of the atmosphere, gradual 
changes occur in the speed of the waves as the temperature, air density, and levels of 
ionization change. On frequencies below 30 MHz, long distance communication is the 
result of refraction of the wave in the ionosphere, where free ions and electrons exist in 
sufficient quantity to affect the velocity of wave travel. Depending on the frequency used 
and the time of day, the ionosphere can support communications from very short ranges 
of less than 100 km (Near Vertical Incidence Signals - NVIS) to distances greater than 
9500 km. 

Ionization of the upper atmosphere is attributed to ultraviolet radiation from 
the sun, and results in several layers of varying densities at various heights surrounding 
the Earth. Each layer has a central region of maximum electron density, which tapers 
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off both above and below that altitude. The ionospheric layers that most influence HP 
communications are the D, E, E s , F„ and F ; layers (Fig. 1). Of these, the D layer, with 
altitudes of 50-90 km, absorbs signals passing through it. The lowest region useful for 
returning radio signals to the Earth at HF is the E layer, with altitudes of 80-150 km. Its 
average height of maximum ionization is about 100 km. The layer refracts waves only 
in the presence of sunlight. Ionization is greatest around local noon, and practically 
disappears after sundown. 




Figure 1 . The Ionospheric Layers [From Ref. 8] 



7 



The sporadic-E region, E s consists of relatively dense patches of ionization that drift 
around from 90-130 km above the Earth Its effects become prominent above 21 MHz 
and into the VHF region 

The region of ionization mainly responsible for long distance communication 
is the F-layer. Its altitudes of ionization range from 150-600 km. It ionizes very rapidly 
at sunrise, reaching peak electron density early in the afternoon at the middle of the 
propagation path. The ionization decays very slowly after sunset, reaching the minimum 
value just before sunrise. During the day, the F region is split into two layers, the F, and 
the F : . The F, layer is usually not an important propagation medium unless it supports 
the only mode to propagate. Its refracting heights are between 150-200 km, forming and 
fading with the passage of the sun. After sunset, the F, layer decays and is replaced by 
a broadened F ; layer, the primary medium supporting HF communications. The thickness 
of the layer ranges from 80 km during the day to a broad 150 km about 320 km above 
the Earth at night. The maximum range of a single hop off the F : layer is about 4000 
km. [Ref. 8] 

Traditionally, the ionospheric layers have been characterized as well-behaved 
stratified layers. From very high resolution measurements of the medium, it is now 
widely accepted that the layers are in continual horizontal and vertical motion Periods 
of high solar activity, with their impact on the ionosphere, often result in unusual 
propagation modes. 
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3. Sporadic E (E s ) 



Sporadic E is a thin reflecting layer in the ionosphere which comes and goes 
sporadically at E-region heights. The most important aspect of E s is the maximum 
electron density or critical frequency, f c E s , and its daily and seasonal variations. [Ref. 8: 
P-29] 

Low-latitude or equatorial E s is basically a daytime phenomenon with little 
seasonal variation. Near the geomagnetic equator the critical frequency f 0 E s exceeds 5 
MHz for 90% of the daylight hours. Equatorial E s is due to a plasma instability caused 
by the high electron drift velocity. It is patchy and transparent, but seems to be a useful 
reflector on long TE circuits, making it possible to get higher MUFs than would be 
expected for normal propagation [Ref. 8: p. 1 25], 

Variations of the ionosphere at low latitudes are so strongly influenced by the 
earth's magnetic field, that it is usually more instructive to consider how the ionosphere 
varies with geomagnetic latitude or with dip angle of the earth's magnetic field rather than 
within a geographic framework. 

4. Trans-equatorial Propagation (TEP) 

In low-latitude ionospheres, the "fountain effect" redistributes electrons, moving 
them from the equator north and south to magnetic latitudes of 10° to 20°. The 
electromagnetic field causes electrons to drift upwards, encountering the horizontal lines 
of force of the earth's magnetic field. Electrons diffuse down these field lines to reenter 
the main body of the ionosphere where the field lines cut through the F region. This 
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causes large clumps of electrons at latitudes 10° to 20° from the magnetic equator. These 
clumps are the peaks or "crests" of the equatorial or Appleton anomaly. [Ref. 8: p.30] 

These crests are most developed in the late afternoon and early evening, during 
the equinoxes, and at solar maximum. Their critical frequencies, f 0 F ; , can exceed 20 
MHz as compared to 10 MHz at the equator. The height of maximum electron density, 
h m F ; , is less at the crests than at the equator. The significant change with latitude of f 0 F, 
and h m F ; is difficult to predict and may be responsible for some of the interesting 
propagation modes observed in the TE environment. [Ref. 8: p.30] 

Discovered in 1947 by radio amateurs, trans-equatorial (TE) propagation occurs 
on circuits which cross the equator, more or less at right angles, and have MUFs 
(maximum usable frequency) higher than the normal multihop modes. Stations attempting 
TE contacts must be nearly equidistant from the geomagnetic equator. Two types of TE 
propagation depend on different features of the equatorial ionosphere for their 
characterization [Refs. 9,10,11,12], These are aftemoon-type TE propagation (A-TEP) 
and evening-type TE propagation (E-TEP). Both types have been observed 
simultaneously on some circuits around 2000 local time [Ref. 8: p. 1 26] . 

a. Afternoon-type TEP 

A-TEP has the following characteristics: 

• MUF greater than the normal 2F MUF, i.e., greater than 40-50 MHz. 

• peak occurrence from 1700-1900 local time, near the equinoxes and at solar 
maximum. 

• path lengths of greater than 6000 km. 
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strong steady signals with low fading rates and small doppler spread 



Several theoretical models have been proposed in the past to describe the 
propagation mechanism in TEP circuits. Initially a double refraction scheme from the 
ionospheric crests was proposed [Ref. 13], Although this model explains quite 
satisfactorily the A-TEP phenomena, it fails to predict the E-TEP basic characteristics 
[Refs. 9, 14]. 

Raytracing algorithms [Ref. 15] have determined the propagation mode 
for these signals is a "super mode" or FF mode, where the signal is reflected twice by the 
F layer, on opposite sides of the equator, without a ground reflection. Figure 2 illustrates 
this super mode, which is dependent upon the electron density concentrations in the 
anomaly crests already discussed. The crests occur at about 15° dip angle north and south 
of the magnetic equator where h m F ; is at a minimum. Therefore, a ray leaving a suitably 
placed transmitter can be reflected from the first crest in a direction to miss the Earth and 
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Figure 2. The "Super-Mode" or FF mode responsible for A-TEP [From Ref. 8] 
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strike the ionosphere at the opposite crest, before reflection back to the receiver The 
MUF would be higher than the normal 2F MUF since the ionosphere at the crests is tilted 
upwards towards the magnetic equator, giving rise to larger angles of incidence. The 
MUFs are also higher because the critical frequencies are so high in the crests. 

Maximum observed frequencies (MOF) for this super mode can exceed 50 
MHz. The high signal strengths are the result of focussing effects as rays arrive from a 
large range of elevation angles Signals also pass through the absorbing D region only 
twice, as opposed to four times for a 2F mode, and are, therefore, less attenuated. This 
supermode occurs most frequently during equinoxes at solar maximum, and then not every 
day. Its occurrence depends on how the crests developed for a specific day. 

b. Evening-type TEP 

E-TEP usually supports higher frequencies than A-TEP and has different 

characteristics: 

• peak occurrence from 2000-2300 local time, near equinoxes and near solar 
maximum. 

• high signal strengths but with deep and rapid fading, and a large doppler 
spread, which can exceed 40 Hz. 

• path lengths shorter than A-TEP, about 3000-6000 km 

• higher MOFs than A-TEP, can be greater than 100 MHz. 

Based on measurements of elevation angles and group delays on a circuit 
between Japan and Australia a waveguide model has been suggested for the E-TEP [Ref. 
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16], The guiding of high frequency waves through field aligned irregularities has been 
examined theoretically by Nielson [Ref. 17], and this theory has been extended by using 
numerical techniques to explain TEP phenomena by Heron and others [Ref. 18], A 
"whispering gallery" mode has also been considered to describe TEP [Ref. 19], Data 
gathered for the fine structure of ionosphere has shown that there are elongated 
irregularities aligned with the geomagnetic field lines in the equatorial ionosphere [Ref. 
20], These are tubular shape depletion regions inside the equatorial ionosphere extending 
on both sides of the magnetic equator at least from 5° to 10° (magnetic dip). [Ref. 14] 
The propagation mode for E-TEP is probably the "whispering gallery" or 
"field-guided" mode. Range spreading on evening ionograms indicate that the equatorial 
ionosphere is threaded with "empty" tubes aligned along magnetic field lines where the 
electron density is much lower than that of the surrounding ionosphere. Figure 3 shows 
that the propagation takes place by rays skidding around the walls of the tube, bouncing 
off the walls, and emerging at the far end of the tube [Ref. 21], 
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MOFs have been noted over 100 MHz, arising from very high angles of 
incidence. Large doppler shifts would be caused by the upward movement of the tubes, 
which rise rapidly to the top of the ionosphere after their creation near the base. The best 
circuits to support this theory place the transmitter tangential to the earth's magnetic field 
so rays commence at the altitude to enter the tubes. Highest MUFs would be achieved 
when the receiver is similarly placed. The circuit, therefore, should be symmetric about 
the magnetic equator, with transmitter and receiver located at magnetic conjugate points, 
as is the case with the Kauai-Rarotonga path. Like A-TEP, E-TEP is unpredictable night 
to night. 
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II. THE AMBIENT COMMUNICATIONS MODEL (AMBCOM) 



A. INTRODUCTION 

This chapter provides a brief description of applicable portions of the AMBCOM 
program. Details of the complete program and code are contained in the user's guides 
published by SRI International [Refs. 22,23,24], 

AMBCOM was designed for batch processing using card images as input. Separate 
programs support each of the function areas modelled, i.e., the ionosphere model is 
generated by a program which passes its data to another program for calculating 
raytracing curves. The system flow is first explained followed by the programs and inputs 
which generated the TE data for the Kauai-Rarotonga path 

B. SYSTEM FLOW 

The AMBCOM system is a multiprogram batch system written in FORTRAN and 
was executed on a VAX 3100 workstation at NPS. Programs NATGEN, RAYTRA, and 
COMEFF were used to generate the TE data, as shown in Figure 4. Data are passed 
between these programs by means of saved data files [Ref. 22: pp. 107-1 55], 

The execution input streams were constructed for each of the 123 days in the period 
of interest, July-October 1962. A file of 57 lines of code comprising the input streams 
for NATGEN, RAYTRA and COMEFF was executed successively for each day. In order 
to obtain a complete 24-hour table of the spectrum of interest, each day's input file 
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Figure 4 . AMBCOM System Flow [From Ref. 25] 



included as a parameter the average of the eight daily Kj, indices (App. B). Total 



execution time for the 123 day campaign was three and half hours. 



C. NATGEN 

1. Overview 

The purpose of NATGEN is to model the ionosphere along the 
communications path between two points. NATGEN builds a model of the F ; , F,, and 
E layers at control points along the path. The control points are evenly spaced at 100 km 
increments with a maximum of 41 control points for paths longer than 4000 km. The 
control points are not necessarily located at the signal reflection points since raytracing 
is performed in RAYTRA. [Ref.22: pp. 23-24] 

NATGEN begins by reading ionospheric coefficients provided by the Institute 
for Telecommunications Sciences (ITS) in Boulder, Colorado, based upon the month, day, 
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and the current sunspot number (input stream variables). AMBCOM uses the ITS Blue 
Deck (a reference to earlier coefficient files that were issued on color coded computer 
cards) as the starting point for modeling the F,, F,, and E layers at each control point. 
Additional parameters describing atmospheric noise, ground conductivity, and ground 
permittivity along the path are read and passed to RAYTRA. [Ref. 22: p.23] 

The semithicknesses, heights of maximum ionization, and the critical 
frequencies for the F 2 , F,, and E layers at each of the control points are passed as outputs 
from NATGEN. ITS upper, median, and lower decile values for sporadic E critical 
frequencies are passed to RAYTRA which performs all E s calculations. NATGEN also 
passes on the location of the transmitter and receiver sites, the time of year, time of day, 
the current K p index, the sunspot number, distance between control points, the number of 
control points, and the path length. [Ref.23: pp. 123-1 24] 

2. Layer modeling 

The ITS ionosphere provides the vertical incident critical frequency (fj for the 
E layer, the F 2 layer, and sporadic E. The E layer critical frequency (4£) is always set 
to the median decile values. For the sporadic E and the F 2 layers, the upper, median, and 
lower decile values represent the critical frequencies 90%, 50%, and 10%, respectively, 
of the days for a given time and month [Ref. 26: p. 89]. NATGEN allows the user to 
choose which of the F 2 layer values will be used as the F 2 critical frequency (4F 2 ) for all 
ionospheric calculations. The default /^ 2 value is the median decile number, and was 
used in this study. Two additional parameters for the F 2 layer are the ratio of the 
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semithickness of the layer (y m ) to the height of maximum density' (h m ) and the maximum 
usable frequency (MUF) for a 3000 km path (M3000). [Ref. 22: p.25) 

Using these coefficients, AMBCOM models the ionosphere at the control 
points. The F,, F, and E layers are represented as three parabolic layers with the values 
for the F, layer derived from the F : and E layer models. The height and semithickness 
of the E layer are set at 1 15 km and 25 km respectively. The F 2 layer height is calculated 
using the ITS coefficients in a two step process. First the peak height of the F ; layer 
(HP f: ) is calculated using the Shimazaki equation [Ref. 22: p.26] 

HP f; = (1490/M3000) -176 (4) 

This value is corrected for signal retardation with a height factor (Ah) equation [Ref. 22: 

p.26] 

Ah = [ f/f c l n { (f/f c + 1 )/(f/f c - 1 ) } -2]y m (5) 

where 

Ah = the height error, 
f c = the critical frequency of the layer, 
f = the transmitted frequency and 
y m = the semithickness of the layer. 

The Ah factor is then subtracted from HP F , to produce a corrected height for the F : layer. 

The F, layer parameters are not represented by ITS data but are calculated 
using the E and F, layer parameters. The F, layer must overlap the F ; layer by half of 
its own semithickness with the bottom of the F, layer set at 130 km. The critical 
frequency of the F, layer is calculated based upon the critical frequency of the E layer 



18 



(f„E). In the event that the F ; layer critical frequency is lower than that for F, . f F is 
set to 0.695. [Ref. 22: p. 27J 

3. NATGEN Input Variables 

AMBCOM is a batch system in which an input stream controls the execution 
of the system and resembles a series of 80 column computer cards. AMBCOM was 
originally designed during the 1970s when card input systems still predominated [Ref. 
22: p.15). Figure 5 is an example of statements used to execute NATGEN. The ASSIGN 
statements provide NATGEN access to the ITS file (ESSABLU.DAT), the conversion 
tables for the geomagnetic coordinate system (RAGCOT.DAT), two output files 
(IONOS.DAT and N01SE.DAT), and execution control cards one through four. 

AMBCOM provides a large amount of flexibility in the control of data 
production. NATGEN models the ionosphere along any propagation path at any hourly 
time increment specified by the numbered control cards. The AMBCOM user defines the 
problem by specifying the transmitter and receiver locations, the sunspot number, the K f 
index, and the time increment. 

Control parameters were chosen to model the ionosphere, as closely as 
possible, to the 1962 SRI transmission paths. The intent was to execute the model from 
the perspective of a communicator who is attempting to estimate the possibility of 
communicating with another HF site under a given set of circumstances. Although the 
communicator does not have the advantage of using the sunspot number and K p index 
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$! 34567 89012345678901234S67B901234S678901234S678901 2345678 901 2345678901234567*90 

SI ••• OCTOBER 1962 ••• 

$1 FREQUENCY MANAGEMENT BETWEEN A SPECIFIED TRANSMITTER AND RECEIVER 
SI INPUT SETUP FOR NATOEN ••• 

SSET DEF DKAO : (AMBCOM ) 

SASSIGN DKAO: (AMBCOMlESSABLU.DAT FOR001: 

DKAO: (AMBCOMlRAGCOT.DAT FOR075: 

DKAO : (AMBCOM) IONOS.DAT FOR002: 

DKAO: (AMBCOM)NOISE. DAT FOR0O4 : 

SYSSINPUT FOR0S8: 

DKAO: ( MCKINSTRYlNATOUT.DAT FOROS9: 

( AMBCOM JNATGEN 

0 0 0 



SASSIGN 
SASSIGN 
SASSIGN 
SASSIGN 
SASSIGN 
SR DKAO 

10 0 0 
BEGIN THESIS WORK WITH AMBCOM 
2 1 

3 22.00 -159 . SO -21.20 

4 0.0 22.0 2.0 



-159.80 



1962 



1 

10.0 



43.0 



4.4 



Figure 5 



NATGEN control cards used to model Kauai-Rarotonga path 



tables for a current month, he should have a rough estimate of the current sunspot number 
and K indices at any given time. 

Figure 5 shows the four control cards used to execute NATGEN. Card one 

is used to edit ionospheric parameters, and card two causes NATGEN to run the 

ionospheric generator program. Card three defines the geographic latitude and longitude 
' \ 

of the transmitter and receiver; negative numbers indicate a west longitude or south 
latitude. Following the location information are the year, month, sunspot number, and K p 
index. Card four indicates the beginning hour, ending hour and time increment. 

Input streams w'ere divided into one-day periods for the four-month study, with 
each day represented by an average K p index. The values for sunspot number and K p 
index were taken from the Journal of Geophysical Research [Refs. 4,5, 6, 7]. 
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I). RAYTRA 



1. Overview 

The raytracing program (RAYTRA) model can be executed in two modes, 
point-to-point and radar. The point-to-point mode, used in this thesis, performs raytracing 
from a transmitter to a receiver. RAYTRA computes group times, phase times, signal 
losses, the effects of sporadic E and elevation angles at both sites. These data are saved 
for each successful propagation path and are passed to COMEFF. 

The raytracing algorithm computes the propagation path based on data 
produced by NATGEN. A raytrace which ends within 1000 km of the receiver site is 
saved for further processing. When two rays bracket the receiver, the program 
interpolates a ray that falls within some preset value, in this case ten kilometers. The 
parameters describing this ray are saved for COMEFF. In the case where only a single 
ray is found, the program again interpolates until a ray close to the receiver is found. The 
AMBCOM User's Guide for Engineers [Ref.22: pp. 35-56] provides an in-depth treatment 
of raytrace algorithms. 

RAYTRA estimates the amount of absorption for each ray. The absorption 
calculations are divided into four parameters as follows: 

• L d is the divergence, i.e., free space spreading loss. 

• L a is ionospheric absorption loss such as deviative, nondeviative, and auroral 

losses. 

• L Es is the loss resulting from sporadic E. 

• L Cj is the loss due to ground reflection along the path. 
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The total loss for a path is the summation of all of these factors. A full explanation of 
the RAYTRA path loss algorithms may be found in [Ref. 22: pp, 56-80], 

2. Sporadic E Calculations 

RAYTRA offers two choices for computing sporadic E (E s ). The first method 
computes the reflection of a signal for frequencies below the blanketing frequency (f b E s ). 
The blanketing frequency is computed based upon the location of the control point The 
second method performs E s reflection calculations for all frequencies less than f 0 E s . The 
value for f 0 E s is based upon the upper, median, or lower ITS f 0 E s coefficients passed from 
NATGEN. The choice of value is specified by the user. RAYTRA calculates f 0 E s 
dependent upon the percentage of E s specified by the user. If 90% E s is specified then 
RAYTRA uses the upper decile value of f 0 E s from the ITS file. For 50% E s RAYTRA 
used the median decile value 

3. RAYTRA Input Variables 

RAYTRA, like NATGEN, provides the user with a number of options for 
controlling program execution. Options are specified with numbered control cards which 
represent input streams following the program execution statement. Figure 6 is an 
example of the RAYTRA controlling statements used for this thesis. 

Card two controls the amount of E s for a given execution. For this study, the 
default value of 50% was used. Card two also allow the user to set the level of man- 
made noise at the receiver. A residential or suburban setting w*as assumed as a liberal 
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$1 34567890123456789012345678901234 5676901234567890123456709012345 67690123 4567090 
$1 * * * INPUT SETUP FOR RAYTRA 

$SET DET DKAO: (AKBCOM] 

$AS5IGN SYS $ I NPUT FOR060: 

$ AS S I GN DKAO : ( MCK INSTR Y J RAYOUT . DAT TOR059: 

$ASS I GN DKA0: ( AMBCOM ] COMEFT . DAT FOR007: 

{ AKBCOM ) RAYTRA 
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1 



22.00 
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Figure 6 . RAYTRA input stream used to model Kauai-Rarotonga path 



estimate of the noise environment. 

Frequencies for RAYTRA are specified by control card three. Frequencies 
from 4-64 Mhz were initially executed to model as closely as possible the SRI 
experiment. Unfortunately, the antenna configuration pattern in AMBCOM only allows 
frequencies of 3-33 Mhz to be transmitted, so the 4-64 Mhz input range always failed 
to show propagation above about 40 Mhz. For the remainder of the study, RAYTRA 
input was set to 4-40 Mhz, with successful propagation up to about 30-32 MHz. 



E. COMEFF 

1. Overview 

The communications effect (COMEFF) program evaluates the quality of a 
transmitted signal. NATGEN and RAYTRA dealt with the ionosphere and possible 
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ionospheric effects along a given propagation path The only site-dependent information 
for those two programs were the site locations. COMEFF parameters include transmit 
power, antenna configuration, signal bandwidth, and other communication site-dependent 
data. Based upon these local parameters and the data produced by RAYTRA, the 
performance of a particular communications link using a given frequency can be 
characterized in terms of the signal-to-noise ratio (SNR), field strengths and doppler 
spread. 

2. Program Features 

COMEFF uses the data generated by RAYTRA for all of the saved modes as 
the basis for evaluation of a particular communications link. The raytrace descriptive 
data, for each ray, includes: 

• take-off angle, 

• group time, 

• phase time, 

• path loss, 

• noise power (per 1 Hz), 

• arrival angle and 

• transmitter frequency [Ref. 22: p. 1 12]. 

Multiple raytraces may be described for each frequency. COMEFF combines the effects 
of multiple modes to produce a single reception statistic. 

COMEFF allows the user to include the effects of the antenna configuration 
on the communications link. The antenna data is taken from ANTLIB.DAT, the 
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AMBCOM antenna file ANT726 was used as a Granger Associates vertically polarized 
log-periodic antenna [Ref. 27], Appendix A contains the antenna gain table for the 
ANT726 antenna used at both transmitter and receiver in this thesis. 

The COMEFF program provides a variety of output options. COMEFF will 
calculate the SNR, group time, doppler shift, phase, and delay spread for all frequencies 
specified in RAYTRA. Other options provide field strengths or analysis of signal quality 
as bit error rate. All requested data may be printed in several formats which include a 
listing of SNR, doppler shift, phase shift, and group times for each mode. COMEFF also 
produces a single report that sums the data for each frequency at a particular time. The 
SNR specified in the COMEFF summary report is a weighted accumulation of SNRs for 
all modes received from RAYTRA. 

Weighted SNRs from the COMEFF summary report were compared to SRl's 
propagation spectrum, which lacked specific SNR data. This composite SNR is calculated 
by computing the power for each ray, combining these values and then subtracting the 
noise value. The power A, for each ray is computed with the equation 

A,= (P l G l G r c 2 )/47tL,f 2 10 i: L r (6) 

where 

P, = transmitter power (W), 

G t = transmitter antenna gain, 

G r = receiver antenna gain, 

L, = path loss for the ith ray, 

L r = receiving antenna loss, 
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c = speed of light and 
f = frequency (MHz). 

The composite SNR value is given by the equation 

SNR = 10 log l0 (ZA,) - N e , (7) 

where N e is the noise power density in dBW. [Ref. 22: p.83] 

3. COMEFF Input Variables 

Figure 7 is an example of the COMEFF input stream. Two control cards 
control program execution. The first card specifies transmit power, and path loss 
threshold. A bandwidth of 4 kHz and a transmitter power of 30 kW for SRI's sounder 
transmissions results in 7.5 W/Hz. The path loss threshold excludes modes with a path 
loss greater than the specified threshold. This parameter was set at 300 dB which 
effectively allowed all modes to be included. The second COMEFF control statement 
specifies antennas. 



$i * * * input setup r OR CONErr 

SI345678901234 56 78901234 56 789012345678901234S678901234 567 89 0123 

$SET DET DKA0 : ( AM8COK ) 

$ ASSIGN SYSSINPUT rOR055: 

SASSIGN DKA0: l MCEINSTRYlHAVrORN.DAT rOR020: 

SASS1GN DKA0: (MCKINSTRYlCOMOUT.DAT rOR059: 

SASSIGN DKAO: [ AKBCOH ) ANTLIB . DAT rOR006: 

SR DRAO : ( AKBCOH ] COHEFT 
RAUA1 -RAROTONGA PATH 25 OCTOBER 1962 

0 1. 2200 010 7.5 0. 300 1 

ANT726 ANT726 1 

SSET DET (HCRINSTRY) 

SEXIT 



45678901234567890 



0 0 1.0 



Figure 7 . COMEFF control input used to model Kauai-Rarotonga path 
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III. THE SRI TRANS-EQUATORIAL EXPERIMENTS 



A. GENERAL 

From June through October 1 962, Stanford Research Institute (SRI) operated oblique 
incidence sounders on a 4800-km TE path from Kauai to Rarotonga as part of the test 
instrumentation for the 1962 nuclear tests in the Pacific. Workers observed anomalous 
propagation of HF modes from low HF into the VHF range across the TE path, and 
noticed this as a nocturnal phenomenon, occurring between sunset and sunrise. The 
propagation was unusual in that frequencies much higher than would usually be predicted 
were propagated over very long distances. This occurred on over 80% of summer nights 
and nearly 100% of nights during the equinoctial months. It was also observed that the 
mode showed a slight inverse correlation with magnetic activity, appearing later on 
magnetically disturbed days. The mechanism involved was believed to depend on field- 
aligned ionization and dip angle as it relates to the magnetic field symmetry. [Ref. 1] 

B. THE EXPERIMENT 

The observations made by SRI resulted as an outgrowth of the 1962 Pacific nuclear 
test series. Specific paths and equipment used were not selected specifically for the study 
of TE propagation. 

One path studied was from Kauai in the Hawaiian Islands to Rarotonga in the Cook 
Islands. The two terminals of this 4800 km path are very nearly geographic and magnetic 
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dip conjugates, with the magnetic equator very near midpath as shown in Fig. 8 The 
great circle path is tilted about 10° from the magnetic meridian. This path will be used 
for the PENEX trans-equatorial experiments. SRI workers used a Granger Associates 
oblique-incidence sounder transmitter at Kauai and a sounder receiver at Rarotonga. The 




Figure 8. Map Showing Location of Kauai-Rarotonga Path [From Ref. 1: p.4] 
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;ounders swept from 4-64 MHz at various time intervals. Peak-pulse power output was 
iO kW, using pulse widths of 100 ps (16 kHz bandwidth) and 1.5 ms (4 kHz bandwidth). 
SRI used vertically polarized log-periodic antennas with 7-8 dB gain over isotropic and 
i beamwidth of 1 10°. [Ref. 1: p.3] 

The sounders operated from June through October 1962 providing data for all but 
i few days of that period. Both the 100 ps and 1.5 ms received signals w-ere recorded 
>n 35-mm film and later digitized to punched card format for their statistical analysis. 
Ref. 1: p.5] 

I. ANALYSIS OF SRI DATA 
1. Time/Frequency Behavior 

The propagation frequency spectrum and the percentage of occurrence of 
>ropagation on any frequency for one-hour periods are shown in Figs. 9-12. Percentage 
iccurrence is indicated by the length of the horizontal line. A line extending completely 
icross a one hour period indicates 100% occurrence during that month. Also shown are 
;unrise and sunset times. [Ref. 1: p.6] 

Figures 9-12 indicate that the VHF mode is a nocturnal process, commencing 
»enerally at sunset, although it occasionally appeared up to two hours before sunset. The 
iverage onset time for July and August was 0500Z, and for October and September, 
1340Z. It was observed that the normal (2F 2 , 3F 2 , M, etc.) modes gradually faded out as 
he VHF mode appeared. [Ref. 1: p.6] 
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Figure 9. Propagation Spectrum, July 1962 [From Ref 1: p 7] 
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Tigure 10. Propagation Spectrum, August 1962 [From Ref. 1: p.8] 
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Figure 11. Propagation Spectrum, September 1962 [From Ref. 1: p.9] 
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The equinoctial month of October (Fig. 12) reveals that propagation occurred 
on almost all frequencies from 5-30 MHz nearly 100% of the time, save for two hours 
after sunrise. Figure 12 also indicates nearly 100% occurrence for 50-64 MHz between 
sunrise and sunset. Seasonally, then, September/October had a higher percentage of 
propagation on most frequencies than did July/August. [Ref. 1: p.6] 

2. Correlation with Spread F and Magnetic Activity 

Geomagnetic activity, as influenced by solar flares and sunspot activity, is 
known to influence ionospheric phenomena such as spread F. Figure 13 shows the 
occurrence of the VHF mode on quiet and disturbed days for the summer months of June, 
July, and August (Fig. 13a) , and the equinoctial months of September, and October (Fig. 
13b). The days were selected from the ten quietest and ten most disturbed days of each 
month based on both and C p , the magnetic character. The curves of VHF mode 




Figure 13 . Occurrence of VHF Mode on Quiet and Disturbed Days during the 
Summer and Equinoctial Months of 1962 [After Ref. 1: p. 14] 



34 



occurrence exhibit many of the same characteristics of spread-F occurrence. There is a 
negative correlation with magnetic activity where onset time is earlier on quiet days than 
on disturbed. Also, seasonally, there are earlier onset times, a higher percentage of 
occurrence, and longer persistence during the equinoctial months than during the summer 
months.[Ref. 1: p.l 3] 

3. SRI Conclusions 

SRI concluded that only two mechanisms for long range TE propagation of 
frequencies greater than about 30 MHz could explain their observed anomalies: super 
mode propagation and the guidance of waves by field-aligned ionization. A super mode 
would rely upon a pair of local maxima of ionization such as presented by the 
geomagnetic or Appleton anomaly. This pair, occurring at ± 20° magnetic latitude, allows 
two ionospheric reflections without a ground reflection This super mode mechanism is 
generally observed between 1200-2000 LT; SRI's TE anomalies were nighttime 
phenomena, occurring at HF down to 4 MHz, so this was probably not the dominant 
mechanism. Field-aligned ionization at F-region heights, however, is an almost nightly 
occurrence. To invoke this mechanism, the magnetic field symmetry seems critical — the 
dip angle must be just right. Since the terminals of the path were within sight of the 
equatorial spread-F belt, good coupling into field-aligned irregularities would be possible 
at altitudes above 200 km with reasonable take-off angles. [Ref. 1: p.33] 
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I). CONVERSION OF SRI DATA 



Data from Figs. 9-12 were converted to surface contour plots for ease of comparison 
with the AMBCOM predicted spectra. Data was taken from each figure at two hour 
intervals since this is the interval for propagation chosen for AMBCOM's COMEFF 
summary report. For each two-hour period, propagation was considered to have occurred 
if the horizontal line which SRI used to indicate percentage of occurrence was greater 
than 507r of the two-hour period. Since amplitude data was not available from Ref. 1. 
contours were constructed for any level of observed propagation. Resulting contour plots, 
representing the conversion data, are shown in Figs. 14-17. This data will be compared 
with the predicted propagation patterns in Chapter IV. 
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Figure 14. Converted Propagation Spectrum, July 1962 
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Figure 15 . Convened Propagation Spectrum, August 1962 
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Figure 16. Converted Propagation Spectrum, September 1962 
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Figure 17. Convened Propagation Spectrum, October 1962 
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IV. AMBCOM DATA ANALYSIS 



A. ANALYSIS METHODOLOGY 

AMBCOM's COMEFF program produces a summary table of composite SNR 
values (in dB), weighted by amplitude of each mode, for all frequencies and times for 
which calculations were desired. Figure 18 is an example SNR summary table for 
October 9, 1962. Frequencies were desired from 4-40 MHz; time (GMT) was plotted at 
two hour intervals. Only those frequencies which supported propagation are presented 
as output. To produce the data necessary for comparison with the contour plots of the 
known SRI data, a selectively sampled set of predicted data was desired. From each of 
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Figure 18. Example SNR table output from COMEFF 
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the four month's SNR data, the five days exhibiting the most magnetically disturbed 
character (D in tables of Appendix B) and the five days exhibiting the quietest magnetic 
character (Q in tables of Appendix B) were grouped separately as samples. SNR values 
of each month's disturbed days were averaged as were each month's SNR values for quiet 
days. 

An appropriate graphical presentation of the mean SNR values was considered to 
be a surface-type contour plot, as other three-dimensional graphs, such as histograms, 
failed to display necessary detail of the data. Data was often not continuous in nature 
and trends were best observed in the contour plot format. 

B. PREDICTED PROPAGATION PATTERNS 

Figures 19-26 are contour plots of the mean values of SNR for the disturbed and 
quiet days for July, August, September, and October 1962. SNR values were grouped in 
ranges of ten dB, the corresponding contours are shaded with the highest values for SNR 
displayed darkest. Comparison of these plots with the converted propagation spectrum 
plots in Chapter three (Figs. 14-17) reveals that AMBCOM poorly predicted observed 
propagation for all months. However, some trends can be discerned that correspond with 
SRI's observed trends. 

1. July Contour Plots 

AMBCOM failed to predict propagation at 4 MHz for all months and only 
rarely above 32 MHz. It was noted (Figs. 19,20) that AMBCOM failed to predict any 
of the observed propagation (Fig. 14) between 0500-1600 GMT. The trends before sunset 
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at 0450 GMT and following sunrise at 1645 GMT are predicted fairly well; however. SRI 
did not observe the predicted propagation at 30-32 MHz around 2400 GMT. 

2. August Contour Plots 

Predicted propagation for August (Figs. 21,22) extended to higher frequencies 
(up to 38 MHz) around 2400 GMT. Again, this was not observed by SRI (Fig. 15) The 
typical null in propagation just before sunrise (1640 GMT) is correctly shown on both 
predicted and actual plots. Figure 21 for the disturbed August days indicates an erratic- 
patch of propagation at 1800 GMT from 20-22 MHz that SRI also observed. The large 
area in Figure 15 showing propagation up to 33 MHz from 0500-1300 GMT is not 
predicted by AMBCOM. 

3. September Contour Plots 

With the equinoctial months, SRI observed (Fig. 16) almost complete 
propagation. AMBCOM. while indicating much more extensive propagation (Figs. 23,24). 
fails to fully predict such strong propagation. Erroneously, AMBCOM predicted 
propagation from 2000-2400 GMT up to 40 MHz, u'hile SRI did not note the same 
behavior. 

4. October Contour Plots 

Again, SRI observed nearly complete propagation from 4-40 MHz for the 
equinoctial month of October (Fig. 17) with a null around sunrise (1620 GMT). 
AMBCOM predicted this gap (Figs. 25,26), but failed again to predict the propagation 
from 30-40 MHz. 
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Figure 19. Predicted Propagation Spectrum, July 1962, Disturbed 
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Figure 20. Predicted Propagation Spectrum, July 1962, Quiet 
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Figure 21. Predicted Propagation Spectrum, August 1962, Disturbed 
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Figure 22. Predicted Propagation Spectrum, August 1962, Quiet 
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Figure 23. Predicted Propagation Spectrum, September 1962, Disturbed 
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Figure 24. Predicted Propagation Spectrum, September 1962, Quiet 
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Figure 25. Predicted Propagation Spectrum, October 1962, Disturbed 
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figure 26. Predicted Propagation Spectrum, October 1962, Quiet 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

A number of conclusions were reached during this study, as listed below: 

• The bibliography provided at the end of the thesis provided a good background 
about TE propagation effects. 

• The SRI sounder TE data were very useful in obtaining an understanding of 
TE propagation modes for a magnetic conjugate path The presence of 
unusual nocturnal modes were shown to exist. 

• This data provided an excellent basis for the comparison of predicted 
MUF/LUF values to observed MOF/LOF data. 

• The SRI data lacked the amplitude information needed for a full evaluation of 
a TE path for PENEX purposes. 

• The comparison of predicted values with observed values indicated significant 
differences, especially during nighttime hours. This would be expected since 
the AMBCOM prediction program, and other prediction programs, do not have 
a routine for the observed TE modes. 

• While ANT726 from AMBCOM's ANTLIB.DAT is a vertically polarized log- 
periodic antenna [Ref. 27], its antenna pattern probably differs from that used 
by SRI, resulting in differences between predicted and observed data. 
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B. RECOMMENDATIONS 



The AMBCOM program is now available from SRI on PC-based software The 
author has not had the opportunity to use this new product. Currently, input control 
parameters must be located in specific columns, and no data entry programs exist to assist 
the user in entering these parameters. If data input were changed from the old control 
:ard approach to an interactive format for data entry, the program would be much easier 
to use. Additionally, graphical output displays should be added. 

Since sizeable differences were found between the predicted and measured results 
for this TE path, it is recommended that prediction programs such as AMBCOM, 
PROPHET, and IONCAP be modified to include the observed TE modes. 
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APPENDIX A. ANTENNA GAIN TABLE 



Table I . GAIN TABLE FOR ANT726, AMBCOM’S VERTICALLY 

POLARIZED LOG-PERIODIC ANTENNA 



OAl* T AS LI fOS A*SCOK'S TE ST I CALL! 90LAAXEED LOO-9IEIODXC 
AJTESSA AST 7 3 4 



PIEQUEVCT. ME 





2 4 


4 


• 


10 


12 


14 


14 


11 


20 


22 


24 


24 


21 


30 


E 






























L 


4.0 -10.0 -10.0 


3 . 4 


3.2 


2.1 


2 . 3 


2.3 


2 . 0 


1.7 


1.3 


1.3 


1.0 


0.1 


0 . 3 


0.4 


1 


4.0 -10.0 -100 


7.3 


7.4 


7 . 1 


4 . 4 


4.0 


3.3 


3.0 


4 . 3 


4.0 


3 . 3 


3.0 


2 . 3 


2.3 


V 


• .0 -10.0 -10.0 


7 . 3 


7.1 


7.9 


7 . 1 


7.7 


7.3 


4.4 


3 . 9 


3.3 


4.4 


3 .9 


3 . 2 


2.9 


A 


10.0 -10 . 0 -10 . 0 


7 .7 


1 . 2 


1.4 


• . 0 


9 . 3 


9 . 1 


0.2 


7 . 4 


4.3 


3.7 


4 . • 


3.9 


3 .3 


T 


12.0 -10 . 0 -10 . 0 


7.4 


1 . 2 


• . 1 


9.3 3 


L0 1 


10.1 


9 . 2 


• .4 


7.3 


4.7 


3 . 1 


4 . 9 


4 . 3 


X 


14.0 -10.0 -10.0 


7 . 1 


7 . 3 


7 . • 


1 . 4 


9 . 4 


9 . 3 


• . • 


• . 1 


7.4 


4.7 


4.0 


3.3 


4 . 9 


0 


14.0 -10.0 -10 . 0 


4.7 


4.3 


4 . • 


7 . 1 


• . 7 


• . 9 


1 . 4 


7 . • 


7.2 


4 . 7 


4 . 1 


3 . 4 


3.3 


■ 


1C .0 -10 0 -10.0 


4 . 3 


3.7 


3 . 9 


4 . 9 


7 . 9 


• . 3 


7 . 9 


7 . 3 


7.1 


4.7 


4.3 


3 . 9 


3.7 




20 .0 -10 .0 -10 0 


4 . 0 


3 . 2 


3 . 2 


4 . 3 


7 . 3 


7 . • 


7 . 3 


7 . 2 


4 .9 


4 . 4 


4 . 3 


4 . 0 


3 . 9 


A 


22.0 -10 .0 -10 .0 


3 . 1 


4.7 


4 . 4 


3 . 7 


4 . 7 


7 . 2 


7.0 


4 . 9 


4 . 7 


4 . 3 


4 . 4 


4 . 2 


4 . 1 


S 


24 . 0 -10 .0 -10.0 


3 3 


4 . 2 


4 . 0 


3 . 0 


4 . 1 


4.7 


4 . 4 


4 . 4 


4 . 3 


4 . 4 


4 . 4 


4 . 3 


4 . 3 


0 


24.0 -10.0 -10.0 


3 . 1 


3 9 


3 . 7 


4 . 7 


3 . • 


4 . 4 


4 . 4 


4 . 4 


4.4 


4 . 4 


4 . 4 


4 . 4 


4 . 4 


L 


240 -10.0 -100 


4 .1 


3 . 3 


3 . 3 


4 . 4 


3.4 


4 . 2 


4 . 2 


4 . 2 


4 . 2 


4 . 3 


4 . 3 


4 . 3 


4 . 3 


E 


30.0 -10.0 -10.0 


4 . 7 


3 . 3 


3 . 0 


4.1 


3.3 


3 . 9 


3.9 


3 . 9 


3 . 9 


4 . 0 


4 . 0 


4 . 0 


4.0 




33.0 -10.0 -10 0 


4 . 2 


2 4 


2 . 0 


2 . 4 


3 .2 


3.4 


3.4 


3 . 3 


3.3 


3 . 3 


3 . 3 


3 . 3 


3.3 


t 


40.0 -10.0 -10.0 


-1.2 


-2.4 


1 . 2 


2.2 


3.1 


3 .4 


3 . 0 


2.7 


2 . 3 


1 .9 


1.4 


1 . 2 


1.0 


a 




























































0 






























E 


40.0 -10.0 -10.0 


-3 . 4 




-3.3 -2.9 -2.3 




-2.3 


-2.2 








• 2 • 1 




0 






























K 


70 .0 -10.0 -10 .0 






-3.4 -3.1 -4.1 


•4*4 


-4 . 3 


-4 . 4 


-4.3 


-4 . 2 


“t « l 




-3.9 


E 






























E 






























t 


•3.0 -10.0 -10.0 


-7.4 


-I . 2 


-1.3 -1.4 -1.3 


- • . 2 




-7 . 7 












































LOSS 99.0 99.0 


0 . 0 


0 . 0 


0 . 0 


0.0 


0 .0 


0.0 


0 .0 


0 . 0 


0.0 


0.0 


0.0 


0 . 0 


0.0 


■»IDTI 340.0 340.0 


110.0 


110.0 


110 .0 110.0 110 .• 


110.0 


110.0 


110.0 


110.0 


110.0 


110.0 


110.0 


110.0 



54 



APPENDIX B. GEOMAGNETIC AND SOLAR DATA TABLES 



The data contained in the following tables is compiled monthly in the Journal of 
Geophysical Research as Geomagnetic and Solar Data, J. Virginia Lincoln, Editor [Refs. 
4-7], Last column annotations denote Q/q • 10 quiet days, Q - 5 quiet days, D • 5 
disturbed days. 

Table II. GEOMAGNETIC AND SOLAR DATA, JULY 1962 





Kp (3-hour intervals) 


Sum 


Ave. 






Day 


1 


2 


3 


4 


5 


6 


7 


8 


Kp 


Kp 


Rz 


Q/D 


1 


2o 


2o 


3- 


3o 


2+ 


2+ 


2- 


2o 


18o 


2.3 


54 




2 


2+ 


lo 


2o 


2o 


1- 


2- 


3- 


3o 


15+ 


1.9 


39 


Q 


3 


3- 


lo 


2- 


2o 


2o 


3- 


2+ 


1 + 


16- 


2 


38 


q 


4 


2o 


4- 


3o 


2+ 


3+ 


2o 


4+ 


5- 


25+ 


3.2 


30 


D 


5 


3o 


4o 


4o 


3o 


2+ 


2+ 


3+ 


2+ 


24+ 


3 


26 


D 


6 


3+ 


3o 


3- 


4- 


2o 


2- 


2+ 


2o 


21- 


2.6 


20 




7 


lo 


i- 


1- 


2+ 


2- 


2+ 


3o 


3+ 


15o 


1.9 


21 




8 


2o 


3o 


3o 


2+ 


2o 


2o 


2+ 


3o 


20- 


2.5 


16 




9 


2+ 


2o 


2o 


1 + 


1- 


1- 


1+ 


2- 


12o 


1.5 


10 


Q 


10 


2o 


lo 


1 + 


0+ 


3o 


3+ 


2o 


2o 


15o 


1.9 


13 




11 


2o 


2- 


3- 


2o 


3- 


2+ 


3- 


2+ 


18+ 


2.3 


19 




12 


lo 


3o 


2+ 


2o 


2o 


2o 


1+ 


1+ 


15o 


1.9 


11 


q . 


13 


3o 


i+ 


2- 


1+ 


3- 


3- 


3o 


3o 


19- 


2.3 


29 




14 


30 


3- 


3- 


2- 


2- 


2+ 


lo 


2o 


17o 


2.1 


33 




15 


3o 


lo 


1- 


1- 


1- 


1 + 


2o 


2o 


no 


1.4 


21 


q 


16 


2o 


1+ 


0+ 


0+ 


lo 


0+ 


0+ 


Oo 


6- 


0.7 


26 


Q 


17 


0+ 


0+ 


0+ 


0+ 


1- 


0+ 


0+ 


lo 


4- 


0.5 


31 


Q 


18 


1+ 


1- 


1- 


1- 


1+ 


1- 


2o 


lo 


8+ 


1 


29 


Q 


19 


1+ 


2+ 


2o 


30 


2+ 


1 + 


3+ 


3- 


18+ 


2.3 


8 




20 


3+ 


4o 


40 


3- 


2+ 


3- 


3- 


2o 


24- 


3 


14 




21 


3o 


3- 


2- 


3- 


2o 


3o 


2+ 


3o 


21 + 


2.7 


23 




22 


3o 


2o 


2+ 


2- 


1- 


1+ 


1+ 


2o 


14+ 


1.8 


23 


q 


23 


lo 


2+ 


3o 


2o 


1 + 


lo 


3+ 


2- 


)6- 


2 


17 




24 


2+ 


lo 


lo 


1- 


2o 


4+ 


5- 


40 


20o 


2.5 


13 




25 


5- 


4+ 


3+ 


2+ 


2- 


1- 


lo 


2+ 


20+ 


2.5 


11 




26 


3- 


6+ 


(r 


6- 


5o 


3+ 


4+ 


4+ 


37+ 


4.7 


9 


D 


27 


So 


4- 


3+ 


3- 


4o 


5+ 


3- 


40 


31- 


3.8 


9 


D 


28 


4- 


4- 


40 


3- 


3o 


4- 


2+ 


3- 


26- 


3.2 


9 


D 


29 


4- 


2- 


3o 


2- 


+ 


3- 


1+ 


2- 


17o 


2.1 


8 




30 


2o 


3- 


2o 


0+ 


i- 


1- 


1- 


1+ 


10+ 


1.3 


7 


Q 


31 


1- 


1- 


1+ 


1 + 


i- 


2o 


4+ 


4+ 


15+ 


1.9 


0 
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Table III. GEOMAGNETIC AND SOLAR DATA, AUGUST 1962 





Kp (3-hour intervals) 


Sum 


Ave. 






Day 


1 


2 


3 


4 


5 1 


6 


7 


8 


Kp 


Kp 


Rz 


Q/D 


1 


5o 


5+ 


5- 


6o 


4- 


3o 


3o 


4o 


35- 


4.3 


0 


D 


2 


3+ 


3o 


lo 


2+ 


2+ 


2+ 


2+ 


2o 


19- 


2.3 


7 




3 


2o 


3- 


3- 


2- 


3- 


3- 


3- 


3- 


20- 


2.5 


7 




4 


lo 


2o 


2+ 


2- 


2- 


2- 


2+ 


1 + 


14o 


1.8 


10 


Q 


5 


2o 


2o 


3- 


3o 


2o 


2- 


1+ 


1- 


15+ 


1.9 


14 


q 


6 


3o 


3+ 


4o 


3o 


4- 


3+ 


3- 


3+ 


26+ 


3.3 


11 




7 


4o 


3o 


2+ 


2+ 


3o 


4- 


5- 


4+ 


27+ 


3.5 


9 




8 


5o 


4o 


So 


4+ 


4o 


4o 


3- 


4+ 


33+ 


4.2 


8 


D 


9 


5- 


3+ 


4o 


3+ 


3o 


2+ 


3- 


3o 


26+ 


3.3 


0 




10 


3o 


3- 


3+ 


2+ 


2+ 


3- 


lo 


lo 


18+ 


2.3 


0 




11 


lo 


1- 


0+ 


0+ 


1- 


0+ 


1- 


Oo 


4o 


0.5 


0 


Q 


12 


0+ 


lo 


1- 


2- 


2- 


0+ 


0+ 


lo 


7o 


0.9 


15 


Q 


13 


l- 


1- 


lo 


1 + 


1- 


0+ 


lo 


3+ 


9o 


1.1 


24 


Q 


14 


2+ 


2+ 


lo 


2- 


2+ 


lo 


3o 


3+ 


17o 


2.1 


40 




15 


3- 


5- 


So 


3+ 


2o 


3- 


3- 


4o 


27o 


3.4 


50 




16 


5o 


3+ 


40 


lo 


1+ 


2- 


4- 


5- 


25- 


3.1 


£0 




17 


4+ 


3+ 


3o 


3o 


3+ 


4- 


4- 


So 


29+ 


3.7 


53 


D 


18 


3+ 


4+ 


5- 


3- 


3+ 


4- 


2+ 


3+ 


28- 


3.5 


43 




19 


4+ 


4+ 


3+ 


3+ 


3- 


2+ 


2o 


2+ 


25- 


3.1 


45 




20 


2- 


i- 


lo 


1 + 


2- 


10 


1- 


2- 


10- 


1.2 


39 


Q 


21 


2o 


i- 


1- 


lo 


lo 


2- 


1+ 


4- 


12o 


1.5 


36 


M 


22 


5+ 


4o 


5+ 


2+ 


3- 


3- 


5- 


5- 


32- 


4 


30 


D 


23 


4- 


3+ 


3o 


3+ 


3+ 


3+ 


3o 


3- 


26- 


3.2 


27 




24 


3+ 


3o 


4+ 


3+ 


4+ 


3+ 


2+ 


4- 


28- 


3.5 


30 




25 


4o 


3- 


3o 


2+ 


2o 


3- 


2o 


3+ 


22o 


2.8 


14 




26 


3- 


2+ 


2+ 


3- 


2+ 


lo 


lo 


lo 


15+ 


1.9 


7 


q 


27 


2+ 


2o 


2+ 


0+ 


Of 


1- 


lo 


1- 


10- 


1.2 


7 


q 


28 


1+ 


2o 


1- 


2- 


2- 


1+ 


0+ 


a+ 


9+ 


1.2 


14 


Q 


29 


lo 


3o 


3+ 


4- 


3- 


3- 


3+ 


3+ 


23o 


2.9 


8 




30 


2+ 


3o 


3- 


2- 


4- 


3- 


4- 


2o 


22- 


2.7 


25 




31 


5- 


5- 


5+ 


4- 


3+ 


1+ 


3- 


2o 


28- 


3.5 


22 


D 
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Table IV. GEOMAGNETIC AND SOLAR DATA, SEPTEMBER 1962 





Kp (3-hour intervals) 


Sum 


Ave. 






Day 


1 


2 


3 


4 


5 


6 


7 


8 


Kp 


Kp 


Rz 


Q/D 


1 


3+ 


2+ 


5- 


4o 


4o 


5o 


5- 


4o 


32o 


4 


48 




2 


3o 


5- 


5- 


5- 


6- 


4o 


4+ 


3+ 


34+ 


4.3 


57 


D 


3 


4o 


5+ 


4- 


5+ 


4o 


4o 




6+ 


38+ 


4.8 


81 


D 


4 


5- 


6- 


4+ 


5- 


4o 


5o 


4- 


3o 


35o 


4.4 


88 


D 


5 


lo 


3+ 


4o 


3+ 


3o 


4- 


3- 


3o 


24o 


3 


90 




6 


4- 


4o 


5o 


4o 


3- 


2o 


3o 


3+ 


28- 


3.5 


82 




7 


4- 


4- 


4o 


3o 


3- 


1 + 


2o 


1+ 


22- 


2.7 


71 




8 


3+ 


4o 


4o 


3o 


2- 


2o 


lo 


2- 


21- 


2.6 


59 




9 


3- 


3o 


i- 


2o 


2+ 


2o 


3o 


3+ 


19o 


2.4 


58 




10 


3o 


4o 


3- 


3o 


2- 


1- 


i- 


lo 


17- 


2.1 


44 




11 


2+ 


i+ 


2+ 


4- 


1- 


1- 


2- 


3- 


15+ 


1.9 


51 


Q 


12 


3- 


3o 


7+ 


5o 


4- 


5- 


7- 


6- 


39- 


4.8 


62 


D 


13 


4o 


3+ 


4- 


3- 


3o 


4+ 


50 


20 


28o 


3.5 


59 




14 


3o 


3o 


4o 


2- 


lo 


2- 


2o 


2+ 


19- 


2.3 


53 


Q 


15 


3o 


4- 


4o 


3- 


i- 


2+ 


3+ 


3o 


23- 


2.8 


60 




16 


2+ 


3- 


i- 


2o 


2o 


2- 


3o 


3+ 


18- 


2.2 


42 


q 


17 


2o 


3- 


2+ 


1+ 


2o 


1 + 


2- 


2o 


15+ 


1.9 


33 


Q 


18 


2- 


1- 


lo 


lo 


lo 


lo 


1 + 


3- 


10+ 


1.3 


35 


Q 


19 


3o 


4+ 


6- 


40 


3+ 


4+ 


6o 


5- 


35+ 


4.4 


24 


D 


20 


5- 


3+ 


2+ 


4- 


2+ 


i- 


lo 


1- 


19- 


2.3 


24 




21 


2o 


3- 


3- 


2- 


3o 


4- 


3- 


3- 


21o 


2.6 


32 




22 


4o 


4- 


3o 


4o 


4- 


4- 


4- 


2+ 


28o 


3.5 


32 




23 


3+ 


4+ 


3o 


3+ 


1 + 


lo 


1 + 


2+ 


20o 


2.5 


41 




24 


2+ 


3o 


lo 


0+ 


1- 


0+ 


2o 


2- 


11+ 


1.4 


54 


Q 


25 


0+ 


1- 


1- 


2- 


2- 


2+ 


2+ 


3- 


13- 


1.6 


53 


Q 


26 


5+ 


6- 


3+ 


2o 


3+ 


5+ 


4o 


2+ 


31 + 


3.9 


55 




27 


3o 


2+ 


2o 


2+ 


3- 


2o 


lo 


2- 


17o 


2.1 


51 


5 . 


28 


3- 


1+ 


1- 


lo 


2+ 


3- 


2- 


3+ 


16- 


2 ^ 


47 


Q 


29 


2o 


3+ 


3o 


3+ 


2+ 


2+ 


3+ 


5+ 


25o 


3.1 


36 




30 


4o 


4+ 


3o 


3- 


3- 


2o 


2- 


1 + 


22- 


2.7 


46 
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Table V. GEOMAGNETIC AND SOLAR DATA, OCTOBER 1962 





Kp (3-hour intervals) 


Sum 


Ave. 




Day 


1 


2 


3 


4 


5 


6 


7 


8 


Kp 


Kp 


Rz 


Q/D 


1 


3+ 


4+ 


5o 


5o 


5+ 


5- 


6+ 


5- 


38o 


4.8 1 


51 


D 


2 


4+ 


4- 


5- 


4- 


4- 


4+ 


3- 


1+ 


28+ 


3.5 


35 




3 


2- 


2- 


4- 


2+ 


3+ 


2- 


2- 


3- 


20- 


2.5 


22 


a 


4 


2+ 


3+ 


3o 


2o 


1- 


2+ 


2o 


3- 


18+ 


2.3 1 


18 


Q 


5 


lo 


2- 


i+ 


lo 


1 + 


3- 


5- 


3+ 


170 


2.1 


13 


Q 


6 


5- 


4- 


4+ 


lo 


3+ 


3- 


0+ 


2+ 


22+ 


2.8 


25 




7 


2+ 


3+ 


i+ 


2o 


1- 


1+ 


3+ 


4+ 


19- 


2.3 


37 


Q 


8 


6+ 


4- 


3+ 


4+ 


3+ 


5- 


4o 


4o 


34- 


4.2 


40 


D 
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